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ABSTRACT 

Shock breakout is the earUest, readily-observable emission from a core-collapse supernova explosion. Ob- 
serving supernova shock breakout may yield information about the nature of the supernova shock prior to ex- 
iting the progenitor and, in turn, about the core-collapse supernova mechanism itself. X-ray Outburst 080109, 
later associated with SN 2008D, is a very well-observed example of shock breakout from a core-collapse su- 
pernova. Despite excellent observational coverage and detailed modeling, fundamental information about the 
shock breakout, such as the radius of breakout and driver of the light curve time scale, is still uncertain. The 
models constructed for explaining the shock breakout emission from SN 2008D all assume spherical symmetry. 
We present a study of the observational characteristics of aspherical shock breakout from stripped-envelope 
core-collapse supernovae surrounded by a wind. We conduct two-dimensional, jet-driven supernova simula- 
tions from stripped-envelope progenitors and calculate the resulting shock breakout X-ray spectra and light 
curves. The X-ray spectra evolve significantly in time as the shocks expand outward and are not well-fit by 
single-temperature and radius black bodies. The time scale of the X-ray burst light curve of the shock breakout 
is related to the shock crossing time of the progenitor, not the much shorter light crossing time that sets the 
light curve time scale in spherical breakouts. This could explain the long shock breakout light curve time scale 
observed for XRO 080109/SN 2008D. We also comment on the distribution of intermediate mass elements in 
asymmetric explosions. 

Subject headings: supernovae: general - supernovae: individual: SN 2008D - hydrodynamics - instabilities - 
shock waves 



1. INTRODUCTION 



The serendipitous discovery of XR O 080109, associated 
with SN 2008D, on 9 January 2008 (|Berger & Soderberg 



20081 [KQng"& Maccarone|[20Q8l [Soderberg et al.||2008 !) has 
allowed us to view a stripped-envelope core-collapse super- 
nova (SN) from its earliest stages, at or near the moment of 
shock breakout from the progenitor star. The radiation burst 
associated with shock breakout is the first electr omagnetic in- 
dicator of a supernova explosion ( |Colgate| 1 968 [ 1 1 974| ) . Shock 
breakout occurs when radiation trapped in the vicinity of the 
supernova shock is able to escape a head of the shock ([Klein 



than that of the typical Wolf-Ray et star, and Soderberg et al. 
argue that this indicates the need for an optically-thick wind 
around the supernova progenitor. Radio observations of SN 
2008D, however, imply a wind mass-loss rate too low for the 
wind to be optically thick at a radius of 7 x 10 ^^ cm ( Soder- 
berg et aL]|2008[ |Chevalier & Fransson||2008| . These con 



& Chevalier|1978t|Ensman & Burrows] 1992t|Matzner & Mc- 



Kee||1999 ). When this happens, the shock transitions from 
a radiation- mediated shock to a hydrodynamic shock ("Katzj 
et al.|2QQ9]). SN 2008D was a normal Type lb supernova 
( [Modjaz et al.| |2009), indicating a compact progenitor lack- 
ing a significant hydrogen envelope. Shock breakout emis- 
sion from such compact progenitors may retain more infor- 
mation about the nature and shape of the supernova driving 
mechanism than breakouts from larger progenitors with intact 
envelopes (see, e.g., [Couch et al.|2()09j ). 

The discovery of X-ray outburst (XRO) 080109 is described 
by [Soderberg et al.| ( |2008| ). The burst lasted about 500 seconds 
and reached a peak Swift XRT count rate of about 7 counts 
s~^ Based on a Comptonized, non-thermal emission model, 
Soderberg et al. conclude that the origin of the XRO is shock 
breakout at a radius of about 7x10^^ cm. This radius is larger 
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elusions w^ere drawn assuming spherical synmietry and raise 
questions about the actual radius of shock breakout. 

The X-ray spectrum of XRO 080109 can be fit reasonably 
well w ith a power law, a b lack body, or a combination of the 
two ( [Modjaz et al.||2009l ). It is thus reasonable to consider 
both thermal and non-thermal sources of emission in a t tempt- 
ing to explain the outburst. [Chevalier & Fransson| ( |2008 1 
posit that a thermal source with a black-body spect rum is 
plausible within the uncertainties of the observations. |Wang 



et al.| p008) argue, however, that bulk-Comptonization will 



scatter thermal photons to higher energies creating a power- 
law spectrum. Although both thermal and non-thermal emis- 
sion sources may be able to explain the shape of the spec- 
trum, neither can account for the timescale of the X-ray light 
curve in a spherically- symmetric geometry. The characteris- 
tic time for both emission types, measured as the full-width 
at half-maximum (FWHM) of the observed light curve is the 
light-crossing time of the progenitor star, R/c, assuming a 
spherically- symmetric shock breakout. The FWHM of the 
XRO, about 100 s ( .Soderberg et al.|2008t[Modjaz et al.|2009l ). 



indicates a spherical breakout radius of ^ 10^^ cm, greater 
than any plausible progenitor radius and well above where the 
progenitor wind could be optically-thick. This is an additional 
contradiction that is difficult to explain with a spherical shock 
breakout model. 

In this work, we describe the characteristics of aspheri- 
cal supernova shock breakout and compare these character- 
istics to the observations of XRO 080109/SN 2008D. Com- 
pounding observational evidence indicates that core-collapse 
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supemovae, especially those involving envelope- stri pped pro- 
genitors (Type Ib/c), are not spherical (see, e.g., [Wang &| 
Wheeler|[2QQ8] ). In the particular case of SN 2008D, polar- 



ization measurements show that the explosion is not spheri 
cal, with dramatic asymmetries in the structure of some line 
forming regions (Maun d et al.||2QQ9[|Goro sabel et al. 2008). 
Theoretically, current models for the explosion mechanism of 
core-collapse S Ne produce inherently aspherical sh ock waves 
dWheeleretal . 2000 , 2002 ; Blondin et al. 2003 ; Burr ows et al. 



2006; Obergaulinger et al.||2006t|Buras et al.|,2 006i Bu rrows 



et al.,2007} . He re we focus on core-co l lapse SNe driven by 
bipolar jets (see [Khokhlov et"aLl[T999l [Couch et al.||2009l). 



as may aris e from a ma gneto-rotational mechanism ( [Wheeler 



[eraL][20QQl [2002: Burr ows et aL][20Q7l ). These models have 
features that may explain many of the observed featu res of 
core-collapse supernovae th at indicate asymmetry ( Khokhlov[ 
let al.[[1999i [Wheeler et al.1|2 000: Wa ng et al.[[20Q l : Hoflich 



et al.[[2001 



[Wang & Whe^eler,|2008^ [Couch etal. [[2009 ) The 



general features of aspherical shock breakout that we discuss, 
however, apply to arbitrarily aspherical shocks, not just those 
produced by bipolar jets. 

The absence of spherical symmetry dramatically modifies 
the observational characteristics of shock breakout and subse- 
quent stages of emission. We assume black-body emission in 
our models and we apply a detector response function appro- 
priate for the Swift XRT and account for X-ray absorption due 
to neutral matter along the line of sight so that we can make 
a direct comparison to the observations of XRO 080109. We 
show that the timescale of the light curve is not set by the light 
crossing time of the progenitor star but by the shock-cro^^mg 
time. This can account for the length of the XRO associated 
with SN 2008D with a Wolf-Rayet star progenitor of reason- 
able parameters. Further, we demonstrate that the spectral 
shape of aspherical shock breakouts is considerably differ- 
ent from that of a single- temperature, spherically- symmetric 
black-body even if thermal emission is assumed. 

Recently Suzuki & Shigeyama (2010) have reported on 
their study of aspherica l supernova shock breakout fr om blue 
supergiant progenitors. [Suzuki & Shigeyama] p010| ) present 
a semi-analytic method for calculating shock breakout light 
curves based on results of 2D hydrodynamic simulations of 
aspherical core-collapse supernovae. They assume, as we do 
in this work, that the breakout emission is thermal and calcu- 
late bolometric breakout light curves. [Suzuki & Shigeyama] 
([2010) find that the asphericity of the explosion and the angle 
from which the explosion is viewed determine the shapes of 
the resulting light curves. We find a similar result in this work. 
Our study, however, is targeted to explaining the observations 
of XRO 080109/SN 2008D and, as such, we calculate X-ray 
spectra and light curves that allow a direct comparison to the 
observations. Also, our simulations are carried out in a more 
compact Wolf-Rayet progenitor star. 

This paper is organized as follows. In Section]2]we describe 
the hydrodynamic simulations of jet-driven SNe. In Section 
]3] we describe our method of modeling the shock breakout 
emission from our simulations. In Section]?] we present the 
resulting spectra and light curves and compare our spectral 
and light curve models with the observations of SN 2008D. 
We discuss our results and give our conclusions in Section J5| 

2. HYDRODYNAMIC SIMULATIONS 

We have carried out high-resolution hydrodynamic simula- 
tions of jet-driven core-collapse supern ovae using the FLASH 
code, version 2.5 ( [Fryxell et al.[[2000] ). We use two progeni- 
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Fig. 1. — Progenitor model density profiles. Model m2rl is shown as the 
solid line and model m7r6 is represented by the dashed line. Model m7r6 is 
simply a stretched version of m2rl; m7r6 is created by increasing the radius 
of model m2rl according to rnew = rorig + 6.17 x 10~^r^^/ . This establishes 
a model with a more extended envelope but leaves the core mass and radius 
practically unchanged. Model m2rl has a mass of 2.5Mq and a radius of 
1.35 X 10^^ cm. Modelm7r6hasamassof 6.8M0andaradiusof 6.5 x 10^^ 
cm. The constant mass-loss rate wind is not shown. 



tor model stars in ou r calcu lations. The first is model slc5a 
from ]Woosley et al.] (11995]). This model is a non-rotating. 



non-magnetic evolved helium star with a pre-supernova ra- 
dius of 1.35 X 10^^ cm (1.93 Rq) and a pre-supernova mass 
of about 2.5 Mq. The second model used is a stretched ver- 
sion of model slc5a. For this model, we stretch the radial 
coordinates at each model grid point of slc5a according to 
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^ + 6.17 X 10 r^;.^ The physical variables, such as 
density and temperature, at each grid point are left unchanged. 
This function then leaves the core mass and radius approxi- 
mately unchanged while extending the envelope of the pro- 
genitor. The mass and radius of this model are 6.8 M© and 
6.5 X 10^^ cm (9.3 Rq), respectively. Both progenitors are sur- 
rounded by a wind with a mass loss rate of 1 .5 x 10~^ Mq yr~^ 
and wind velocity of 1000 km s~^ . The wind is assumed to be 
spherically symmetric, which may not be the case for a ro- 
tating progenitor. The transition from the progenitor model 
profile to the wind profile is made linearly over about eight 
computational zones. Figure [T] shows the density profiles of 
the two progenitor models. 

We use an equation of state (EoS) that accounts for con- 
tributions to the internal energy and pressure from radiation 
and gas. In the wind, the gas and radiation are not in thermal 
equilibrium and so radiation will not contribute to the pres- 
sure or internal energy, however our single-temperature code 
cannot correctly account for this using an EoS that calculates 
contributions from both radiation and gas. Therefore, we ini- 
tially set the temperature in the wind to a small value. This is 
justified because for a strong shock the upstream temperature 
is unimportant to the downstream thermodynamics. We track 
seven atomic species in our simulations: ^He, ^^C, ^^O, ^^Ne, 
^^Mg, ^^Si, and ^^Ni. We employ the consistent multi-fluid 
advection scheme implemented in FLASH ( ]Plewa & Muller 
]1999] ). We do not include nuclear burning. Gravity is calcu- 
lated using the multipole Poisson solver with m = and l=\. 

All simulations are carried out in 2D cylindrical geometry. 
In order to cover the enormous dynamic range from the inner 
regions of the progenitor star (~ 10^ cm) to the shock radius 
several minutes after shock breakout (~ 10^^ cm), we have 
implemented a logarithmically-spaced cylindrical mesh. This 
is achieved through a radially-dependent maximum level of 
refinement limiter. This limiter requires that the grid spacing 
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at radius r, Ax, not fall below some fraction of the radius r. 
The grid spacing then takes the form Ax > r]N~^r, where A^^ 
is the number of zones per block in the x-direction and r] is 
a small number that sets the resolution scale. In effect, r]N~^ 
is analogous to the minimum angular resolution in spherical 
geometry. Additionally, we have set the maximum level of re- 
finement anywhere on the grid to be time-dependent; succes- 
sively higher levels of refinement are dropped from the grid 
as the simulation proceeds. This has the effect of dramati- 
cally increasing the Courant-limited time step at late times, 
allowing the calculations to be completed in a relatively small 
amount of computer time. Each simulation described in this 
paper required approximately 3000 CPU hours to cover 10^ 
seconds of simulation time. This also negated the need to 
re-map the simulation ont o a new grid to continue the simu- 
lations to la te times (e.g., [Couch et al.|2009; Kifonidis et al. 
[20031 [2006] ). 

The jets that drive the explosions are introduced as time- 
dependent boundary conditions at the inner boundary of the 
grid where we inject two identical, oppositely-directed ener- 
getic flows. In order to facilitate this, an essentially spherical 
inner hole is excised from the 2D cylindrical grid. Within this 
hole, the hydrodynamic solution is not calculated. A diode 
boundary condition was enforced at the edge of the hole (see, 
e.g., Zing ale et al.|2002 ). This boundary condition is equiva- 
lent to an outflow boundary condition when the flux into the 
hole is positive, but the flux out of the hole is always zero. We 
include the gravitational effect of the mass initially residing 
within the hole as a Newtonian point-mass at the center of the 
grid, and compute the self-gravity of the gas on the grid. The 
mass that flows into the hole is tracked and included in the 
calculation of the central point-mass gravitational potential. 
The radius of the hole expands during the simulation, cutting 
out the smallest zones where the Courant condition is most 
limiting and ensuring that the hole radius is always resolved 
by a large number of zones as the maximum allowed refine- 
ment level is reduced. The jet injection velocity, Vjet, varies in 
time according to 



TABLE 1 
Simulation parameters 



t < 0.25 ^jet, 
0.25 /^jet < t < ^jet, 



v^jet(0 = { Vniax^(l -^Ajet), 0.25 /^jet < t < ^jet, (1) 

^0, t > /jet, 

where v^ax is the maximum jet injection velocity and /jet is the 
total jet injection time. 

We ran a total of six simulations. Two of these simulations 
are spherical, non-jet-driven, explosions for comparison to the 
jet-driven cases. The spherical explosions are initiated in an 
identical manner to the jet-driven cases: injection of energetic 
material, except that the "jet" opening half-angle is 7r/2. For 
the four jet simulations, the opening half-angle of the jets is 
about 7r/12. The parameters of the jets are listed in Table [T] 
The model name labeling scheme is mMr7?[cold, hot], where 
M is the progenitor mass to the nearest solar mass, R is the 
progenitor radius in units of 10^^ cm, and the cold or hot des- 
ignates the jet parameters used, given in Table [T] For all sim- 
ulations, the maximum extent of the grid is 10^ cm and the 
initial radius of the inner hole is 2 x 10^ cm, roughly the radius 
of the iron core of the progenitor models used. The ambient 
density and temperature at this inner radius for both progen- 
itors is about 5.2 x 10^ g cm~^ and 3.3 x 10^ K. The jets are 
assumed to consist entirely of ^^Ni to facilitate the tracking of 
the injected jet material. We note, however, the jet parameters 
in some of our models, e.g., m7r6cold and m7r6hot, would 
predominantly freeze out into lighter nuclei (e.g., ^He) and 



Model 


Vmax 


Pjet ^ 


Tjet ' 


H.' 


Mtot ' 


Etot ' 


m2rlcold 


3.3 


25.0 


3.0 


2.00 


0.10 


0.8 


m2rlhot 


1.0 


70.0 


8.0 


2.00 


0.12 


1.4 


m2rlsph 


1.0 


70.0 


8.0 


0.08 


0.09 


1.4 


m7r6cold 


7.0 


3.0 


3.0 


8.00 


0.10 


3.8 


m7r6hot 


2.1 


7.0 


6.2 


8.00 


0.07 


3.7 


m7r6sph 


2.1 


7.0 


6.2 


0.32 


0.07 


3.7 



^ Maximum injection velocity of jets in units of 10^ cm s~^ . 

^ Density of the injected material in units of 10^ g cm~^. 

^ Temperature of the injected material in units of 10^ K. 

^ Total injection time in seconds. 

^ Total mass injected in solar masses. 

^ Total injected energy in units of 10^^ erg. 

not into iron group elements (see, e.g., Pruet et al. 2004). 
Our slower, denser jets would freeze out into the iron group. 
The true resulting ^^Ni fraction will then be a strong func- 
tion of the proton fraction Y^ in the jet that we do not attempt 
to model. The parameters of the simulations were chosen so 
that in every case, the injected jet mass is about O.IM©, a 
value simil ar to the ^^Ni mass estimated from observations of 
SN 2008p fSoderberg et al.|2008tpazzali et al.|2008lpod 



jaz et al.||2009| ). Also, for each model, except m2rlhot and 
m2rlsph, the ratio of explosion energy to ejecta mass is about 
0.8 (10^^ erg/M©), similar to the ratio estimated from mea- 
surements of the photospheric velocity of SN 2008D at maxi- 
mum Hght fSoderberg et al.'2008';'Mazz aH et al.|20Q8| . Model 
m2rlhot and m2rlsph have slightly higher explosion energy 
to ejecta mass ratios of about 1.4 (10^^ erg/M©). The jet 
parameters for models m2rlcold and m2rlhot approximately 
correspond to the jet parameters used in Couch et al. (2009 ) 
for their models v3ml2 and vlml2, respectively. There are 
25 levels of refinement at the start of each simulation and the 
effective angular resolution, r]N~^, is 7r/1024. The simula- 
tions are run until 10^ seconds, long after shock breakout in 
all cases. 

Figures [2] -[5] show show density plots of the four jet-driven 
explosion simulations at three epochs: when jet injection 
stops, initial shock breakout, and the end of the simulation. 
In each jet explosion simulation, the jets drive bipolar shocks 
that expand out from the jet injection sites along the cylindri- 
cal axis. The shocks cross in the equatorial plane and establish 
a dense, hot pancake of unbound material. The shocks in all 
cases erupt from the surface of the progenitor stars first at the 
poles. The shocks accelerate into the low-density wind re- 
gion and sweep around the surface of the progenitor and cross 
again on the equatorial plane. This happens just before the 
original equatorial shock structure erupts from the progenitor 
surface. The prolate shock structure evolves toward sphericity 
in the wind region as the reverse shock, established by the out- 
going shock colliding with the wind, sweeps up an unstable 
shell of ejecta. 

The explosions in the smaller progenitor, models m2rlhot 
and m2rlcold, reach the surface of the progenitor approxi- 
mately 50 seconds after the start of the simulations. For ex- 
plosion m2rlcold, the shocks take about 30 seconds to cross 
the surface of the progenitor and collide along the equatorial 
plane. The shock surface-crossing time is only 20 seconds in 
model m2rlhot because the shock structure is more spheri- 
cal than in m2rlcold. The shocks reach peak speeds of about 
1.4 X 10^^ cm s~^ immediately following eruption from the 
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Fig. 2. — Density plots for model m2rlcold at 2, 46, and 10^ seconds, from left to right. The left panel shows m2rlcold at the time jet injection ceases. The 
bipolar shocks are beginning to cross in the equator, establishing a hot, dense, outward-moving pancake of ejecta. The edges of the jets present Kelvin-Helmholtz 
ripples while the contact discontinuity at the jet head is beginning to show growth of Rayleigh-Taylor fingers. The middle panel shows the simulation at the 
moment the polar shocks are erupting from the progenitor's surface. As the shocks erupt from the surface, they accelerate to speeds approaching 0.5c and quickly 
sweep across the circumference of the star crossing again in the equatorial plane. The right panel shows m2rlcold at the end of the simulation. A thin, unstable 
shell has formed at the contact between the wind and ejecta and has broken up into several Rayleigh-Taylor fingers. The shock structure is slightly prolate with 
an axis ratio of about 1.2. 
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Fig. 3. — Density plots for model m2rlhot at 2, 50, and 10^ seconds, from left to right. The left and middle panels clearly show the greater amount of lateral 
expansion of the jets as compared with m2rlcold (Fig. [2J. The left panel also shows the dense pancake of ejecta formed in the equator by the crossing bipolar 
shocks, as well as the very unstable nature of the contactsoetween the jet flow and the star. The middle panel shows m2rlhot at the first instance the shocks erupt 
from the progenitor surface. The shocks sweep across the progenitor circumference at speeds nearing 0.5c and cross again in equator. The right panel shows 
m2rlhot at the end of the simulation. The contact surface between the wind and the ejecta shows significant growth of RT fingers. The shock structure at the end 
of the simulation is close to spherical. 



progenitor surface and then begin to slow in the wind. At 
the end of the simulations, around one day after shock erup- 
tion, the pole to equator axis ratio for m2rlcold is 1.2 and for 
m2rlhot is 1.0. 

The explosions in the larger progenitor take significantly 
different amounts of time to reach the surface. Model 
m7r6cold takes about 225 seconds to erupt from the progen- 
itor poles, while model m7r6hot takes about 390 seconds to 
do the same. The time it takes the shocks to sweep across 
the progenitor surface is also different, taking 225 seconds in 
model m7r6cold and only 125 seconds in model m7r6hot. As 
is the case for the smaller progenitor simulations, this is be- 



cause the hot-jet model is less prolate than the cold. In the 
larger progenitor, the shocks reach peak speeds after breakout 
of about 1 X 10^^ cm s~^ 

We have carried out a resolution study using jet explosion 
model m2rlcold. We have run two additional simulations at 
resolutions of r]N~^ = 7r/788 and r]N~^ = 7r/1331. The results, 
compared with those of the fiducial simulation with r]N~^ = 
7r/1024, are shown in Figure|6] The higher resolution simula- 
tion shows large-scale of north-south asymmetry. This is due 
to small-scale north- south asymmetries near the cylindrical 
axis early-on that then propagate to large scales as the simula- 
tion proceeds to later times. Artificially accelerated growth of 
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Fig. 4. — Density plots for model m7r6cold at 8, 228, and 10^ seconds. The left panel shows m7r6cold at the time jet injection has ceased. The bipolar shocks 
have crossed in the equator creating a dense pancake of ejecta. The lower-density jet material has already been largely shredded by instabilities. The middle panel 
shows m7r6cold at the time of initial shock eruption. Instabilities in the jet material have created high-velocity fingers that have impinged upon and distorted the 
shock structure. Upon shock eruption, the shocks sweep out into the wind at speeds around 0.3c and cross again in the equatorial plane. The right panel shows 
m7r6cold at the end of the simulation. As in the simulations in the smaller progenitor, a dense, unstable shell has formed at the contact between the wind and 
ejecta. The overall shock structure is very slightly prolate. 
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Fig. 5. — Density plots for model m7r6hot at 8 seconds, 390 s, and 10^ s, from left to right. The left panel shows model m7r6hot at the moment jet injection 
is shut off. Highly unstable and turbulent flow is already evident. The bipolar shocks have crossed in the equator creating a dense toroidal outflow there. The 
middle panel shows this model at the moment of polar shock breakout. This model is significantly less prolate than its counterpart, m7r6cold (Fig.|4j, leading to 
an equatorial shock breakout that is relatively close in time to the polar shock breakout. The right panel shows the density field at the end of the simulation. The 
final structure is slightly oblate. A thin, unstable shell is evident at the contact discontinuity between the supernova ejecta and swept-up wind material. 



instabilities near the axis is a well-known problem in Eulerian 
calculations carried-out in curvilinear coordinate systems and 
higher-resolution simulations are more susceptible to these ar- 
tificial instabilities (see, e.g. Fryx elletal.|1991j ). Additionally, 
it has been documented that the directionally- split piecewise 
parabolic method for Eulerian hydrodynamics does not con- 
serve symmetries in small-scale struc tures (see, e.g., |Liska^ 
|Wendroff|2Q03]|Almgren et al.|2006| ). The higher resolution 
simulation produces slightly higher shock velocities momen- 
tarily during shock breakout, reaching a maximum speed of 
1.5 X 10^^ cm s~^ Within the accuracy of velocity measure- 
ments, the low-resolution simulation attains breakout shock 



speeds equivalent to those of the medium-resolution simula- 
tions, about 1.4 X 10^^ cm s~^ . The shock velocities in the dif- 
ferent resolution simulations quickly become equivalent, as is 
evident by the similar extent of the shock structures shown in 
Figure [6] Because of the influence of the numerical artifacts 
that appear in the high-resolution simulation, the results found 
in the low- and intermediate-resolution simulations are more 
reliable. In fact, in terms of calculating the X-ray emission 
from the simulations, there is negligible difference between 
the low- and intermediate simulations, as shown in Figure [19] 
indicating that the gross dynamics have effectively converged 
at the intermediate resolution. The small-scale differences be- 
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tween the low- and intermediate-resolution simulations have 
little impact on the resulting X-ray emission. 

3. SHOCK BREAKOUT EMISSION CALCULATION 

In order to understand the observational effects that aspher- 
ical shock breakout would produce, we simulate the X-ray 
emission from our simulations. We calculate the X-ray spec- 
tra at each output time from our simulations and produce 
X-ray light curves and integrated spectra. We also compute 
the estimated X-ray counts as would be detected by the Swift 
XRT, accounting for the XRT detector response function and 
photon absorption along the line of sight. In this section we 
discuss the details of our approach. In Section |4] we present 
the simulated spectra and light curves for the explosion sim- 
ulations and discuss the general observable characteristics of 
aspherical shock breakout. 

We calculate the spectra of our simulation results in the 
Swift XRT bandpass (0.1 to 10 keV) as a post-processing step. 
The AMR data are first merged onto a uniformly- spaced 2D 
cylindrical grid using volume- weighted averaging. In order 
to calculate the emission as seen from multiple viewing an- 
gles, the data are rotated prior to calculating the optical depths 
(see Figure [7]). The electron scattering and absorption optical 
depths, Tes and Tabs, are then calculated via integration along 
rays directed from the assumed location of the observer. We 
assume the observer's line of sight to be along the cylindri- 
cal /^-coordinate in the post-rotated data. This line of sight 
is then not normal to the simulation data symmetry axis for 
non-zero viewing angles, as shown in Figure |7] We define 
the thermalization depth, where the radiation and matter tem- 
peratures equilibrate, to be where the effective opti cal depth 



n = \/3Tabsrtot = 2/3, and rtpt = Tes + Tabs (see e.g., |Rybicki 
& Lightman 19861 Ensman & Burrows 1992). We assume 
a black-body emission spectrum is formed with temperature 
equal to the matter temperature at the thermalization depth. 
The emission directed toward the observer is then the black 
body flux from the thermalization depth times the surface area 
of the thermalization depth projected toward the observer. 

The thermalization depth, and hence the emissivity, is 
strongly dependent on the scattering and absorptive opaci- 
ties. We use multi-group opacities obtained from th e TOPS 
database m aintained by Los Alamos National Lab ( [Magee 
|et al.|[T995j ). These opacities are temperature- and density- 
dependent and use 32 photon energy groups spaced logarith- 
mically in the XRT bandpass. We assume a chemical mixture 
comprised mostly of helium, as would be relevant to the Type 
lb SN 2008D, with a solar mix of metals at half of the so- 
lar metal abundance (i.e., X = 0, F = 0.992, and Z = 0.008). 
The inclusion of metals is critically important to obtaining 
accurate values for the absorptive opacities, as we discuss in 
Section |4.1[ The TOPS database gives values for the num- 



ber of free electrons per atom, A/^, the Rosseland mean opac- 
ity, and the Planck mean opacity at each temperature, density, 
and photon energy. We assume the total absorptive opacity, 
^abs(^7). to be the Planck opacity. The electron scattering 
opacity is n^^ = O.IA/^ cm^ g~^ for a predominately helium 
gas. The optical depths are then 



nes,abs] (^7 ^ ^) = / ^[es,abs] (^7 ' ^' ' ^)p(^' ' ^) ^^' ' 
J 00 



(2) 



where R and z are the cylindrical radius and height, and E^ is 
photon energy. A schematic diagram of the optical depth inte- 
gration is shown in Figure [7] The radius of the thermalization 



depth as a function of photon energy and z is then 

/?*(^^,z) = 7?(where r*(£^,z) = 2/3), (3) 

and the temperature at the thermalization depth is T^{E^^z) = 
T(R.). 

The intensity of the X-ray emission is assumed to have a 
black-body spectral energy distribution with a color tempera- 
ture equal to T^(E^^z), i.e.. 



I(E^.z)-- 



24 



1 



y/kT^ - 1 • 



(4) 



The emergent flux is then the intensity multiplied by the area 
projected toward the observer. The projected area in cylindri- 
cal geometry for a line of sight perpendicular to the symmetry 
axis is just Ap(E^^z) = 2R^^o(E^^z)dz, where R^o is the trans- 
formation of R^ into the frame in which the simulation sym- 
metry axis corresponds with the cylindrical axis. For other 
viewing angles the projected area becomes 

Ap(E^,z) = 7rsm^a \R^^o\dR + 2fcos^ a R^^odz, (5) 

where R^o is always measured from the symmetry axis (see 
FigurelT]), / is for R^q < and 1 otherwise, and a is the an- 
gle between the symmetry axis and the z-axis, or equivalently 
the angle between the line of sight and the normal to the sym- 
metry axis. In this way, a surface area element is treated as 
a cylindrical ring and we account for shadowing effects. We 
have assumed that the effective optical depth, r* , is axisym- 
metric. This is only approximately correct as it neglects limb- 
darkening effects at latitudes away from the symmetry axis. 
The right panel of Figure [7] shows a graphical representation 
of the integration for a non-zero viewing angle a. 

The total specific luminosity projected toward the observer 
is 

Lp(E^) = ^I(E^,z)Ap(E^,zy (6) 

z 

The total luminosity directed toward the observer is 
J Lp(E^) dE^. In the calculation of Lp(E^), we correct for 
light- travel time effects. For the arbitrary geometries we con- 
sider, this is accomplished by assuming that the observed 
time cadence is the same as the simulation output cadence, 
^obs _ ^sim _ ^^ Each emitting surface area element is assumed 
to have a constant luminosity in the interval dtn. The arrival 
time of the energy emitted by a surface area element in a given 
time interval, i.e., Lp(E^)dtn, is calculated based on the ge- 
ometry of that area element, and the amount of energy the 
observer would see in time dtn is summed-up. The observer's 
measured luminosity at time tn is then this energy divided by 
the time interval dtn. This approach to correcting for light 
travel time is applicable to arbitrary geometries of emitting 
surfaces. 

The emitting regions of our simulations are typically sam- 
pled by about 1000 lines of sight along which the specific 
luminosities are calculated. We restrict our analysis to the 
Swift XRT bandpass, 0.1 keV to 10 keV. Since the absorptive 
opacity is photon energy dependent, the thermalization depth 
is different for each photon energy group. The total spectrum 
is thus not a single temperature black body but the superposi- 
tion of many black bodies at different temperatures and with 
different emission areas. In order to make a more direct com- 
parison with the XRT observations, we convolve our model 
X-ray spectra with the XRT detector response function and 
account for X-ray absorption due to neutral matter along the 
line of sight. We assume a distance to SN 2008D of 31 Mpc. 
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Fig. 6. — Density plots for model m2rlcold at three different resolutions at 1000 seconds. The extent of the shock structures are very similar for all three 
resolutions. The amount of North-South asymmetry is greater in the high-resolution simulation, due in large part to the influence of the symmetry axis (see text). 
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Fig. 7. — Schematic diagram of the integration to find the optical depths. The gray-scale gradient represents the density and opacity fields. The integration 
of the optical depths is carried out along lines-of-sight from the observer, indicated by the top arrows in the diagram. The pair of dashed lines demonstrate 
the photon energy dependence of the effective optical depth, r* . The right panel shows the integration geometry for a viewing angle not perpendicular to the 
simulation symmetry axis. The shapes of the emitting surface will be different in the rotated case than in the non-rotated case. The right panel also demonstrates 
how R^i^ can be negative. 



4. SIMULATED SPECTRA AND LIGHT CURVES 

We have calculated shock breakout X-ray spectra and light 
curves at various observer viewing angles for the four jet- 
driven explosion models. The light curves are presented in 
Figures 12] - [15] and the time- averaged spectra in Figures [17 
- [18] All emission models are calculated assuming a neutral 
matter column depth along the line of sight, A^h, of 1.7 x 10^^ 
cm~^, corre sponding to the G alactic value in the direction of 
SN 2008D ( [Dickey & Lockm an 1990). X-ray photon absorp- 
tion due to neutral matter could occur both in the host galaxy 
(but well beyond the scales we simulate) and locally in the 
Milky Way. All models are calculated using opacities appro- 
priate for a helium gas with a solar mix of metals with metal 
abundances that are half of the solar values (i.e., Z- O.SZ© 
= 0.008). The influence of varying both A/h and Z are dis- 



cussed in Section |4rT] The general light curve characteristics 
are given in Table[2] 

We show the XRT data of XRO 080109 in Figures [8] and 
[9] To prepare these data, we downloaded Swift observation 
00031081002 from the High Energy Astrophysics Science 
Archive Research Center ^ and reprocessed the XRT data with 
the "xrtpipeline" task using the latest calibration files. We ex- 
tracted events from SN 2008D using xselect and produced a 
spectrum and light curve. Because SN 2008D was piled up 
during this observation, we used an annular extraction region 
of I'.'^l inner radius and 40^^81 outer radius, corresponding 
to the 40% encircled energy radius of the PSF and the 90% 
encircled energy radius, respectively. The response files were 
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TABLE 2 
Light Curve Characteristics 



Model 


Angle ^ 


^X,max ^ 


FWHM'^ 


Afd 


Energy ^ 


m2rlsph 




12.7 


7.5 


137.3 


0.52 


m2rlcold 





3.6 


103.6 


581.4 


0.90 


m2rlcold 


7r/4 


3.7 


123.2 


584.8 


0.91 


m2rlcold 


^/2 


3.2 


117.8 


658.4 


0.85 


m2rlhot 





7.4 


5.9 


518.4 


1.32 


m2rlhot 


7r/4 


7.5 


111.6 


476.8 


1.37 


m2rlhot 


^/2 


4.9 


96.8 


559.2 


1.41 


m7r6sph 




9.4 


27.9 


245.1 


0.64 


m7r6cold 





1.6 


223.5 


825.4 


0.60 


m7r6cold 


7r/4 


1.3 


280.7 


924.0 


0.53 


m7r6cold 


^/2 


0.65 


16.9 


1104.3 


0.34 


m7r6hot 





3.8 


154.9 


235.6 


0.56 


m7r6hot 


7r/4 


1.9 


199.8 


623.2 


0.49 


m7r6hot 


^/2 


3.0 


7.5 


383.8 


0.29 



^ Viewing angles of represents a line of sight along the equator and 

viewing angles of 7r/2 represents a line of sight along the axis of 

symmetry. 

^ Maximum X-ray luminosity in units of 10^^ ergs s~^ . 

^ Full width at half maximum X-ray count rate of the light curve in 

seconds. 

^ Total time over which the X-ray count rate was greater than 10% 

of maximum. 

^ Total radiated energy integrated over At in units of 10^^ ergs. 
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Fig. 8.— The XRT light curve of XRO 080109. The light curve is con- 
structed using a photon binning length of 15 seconds. The FWHM of the 
light curve is about 100 seconds and the rise time to maximum is about 70 
seconds. The X-ray burst lasts about 600 seconds. 




Fig. 9.— XRT spectrum of XRO 080109 integrated over the total burst 
time of 600 seconds. 



produced using the "xrtmkarf " task, which took into account 
the fact that we extracted events from only 50% of the PSF. 
These response files were also used in the post-processing of 
our simulation results to compare the models directly to the 
2008D data. Therefore, the simulation light curves and spec- 
tra can be thought of as having been observed with only ~l/2 
of the XRT effective area. 



Figures 12-15 show both the total X-ray luminosity and the 



predicted X-ray count rate, corrected for detector response. 
X-ray absorption, and distance. As can be seen, the shapes of 
the count rate curve and the luminosity curve are not the same. 
This is because the underlying spectrum is varying in time. At 
times when the spectra are softer, more photons are produced 
at lower energies where the effects of X-ray absorption and 
detector response are stronger. Thus it is not appropriate to 
use a constant count rate-to-luminosity conversion factor at 
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Fig. 10. — XRT band (0. 1 - 10 keV) light curves for the spherical explosion 
models, m2rlsph (solid blue line) and m7r6sph (dashed red line). The time 
coordinates have been adjusted so that the peak of both light curves lies at 
t = 0. The top panel shows the X-ray luminosity and the bottom panel shows 
the XRT count rates corrected for detector response and absorption along the 
line of sight. The width of the light curves is set by the light crossing time 
of the progenitors. The light curve of m2rlsph is comprised of significantly 
harder photons than that of m7r6sph which is the reason the count rate light 
curves show more disparity than the luminosity curves. 

all times. 

In each explosion model, the spectrum varies significantly 
throughout the burst. This is true even in the spherical ex- 



plosions (see Figure [TO]), though to a lesser extent, because 
as the radius of the thermalization depth increases and the 
temperature there cools adiabatically, the reduction in emis- 
sivity is somewhat balanced by the increased emitting surface 
area. Figure [16] demonstrates the temporal variability of the 
spectrum for model m2rlcold. This figure shows the instan- 
taneous spectrum of m2rlcold at four times during the first 
150 seconds of breakout emission. The thermalization layer 
in each explosion lies in the region between the forward and 
reverse shocks. At later parts of the X-ray bursts, the thermal- 
ization layer coincides roughly with the contact discontinuity 
between the swept-up wind and the reverse- shocked ejecta, 
where the density increases significantly. Therefore, the na- 
ture of the breakout emission can depend significantly on the 
character of the wind. The absence of a progenitor wind (or 
the limit of a very low-density wind) may result in the ther- 
malization layers being driven deeper into the ejecta, possi- 
bly below the contact discontinuity. Since the temperature is 
highly discontinuous at the contact discontinuity, this could 
dramatically change the emission. 

The light curves for these jet-driven explosions may be gen- 
erally described as having two peaks. The first is associated 
with the initial shock breakout at the poles. The dip in the 
light curves following this first peak is due to the decreasing 
temperature of the thermalization layer as the shock moves 
out further into the wind. The second peak is attributed to 
emission from the equator as the bipolar shocks cross there 
and create a pancake of very hot, twice-shocked gas. As 
the explosions evolve more toward sphericity at later times, 
the light curves decay exponentially. The time scales of the 
light curves are related to the shock crossing times of the pro- 
genitors. Figure [TT] shows the post-breakout evolution of the 
temperature and thermalization depth at 0.75 keV for model 
m2rlcold. Before shock breakout the location of the thermal- 
ization layer is the transition region between the wind and the 
progenitor star, where the density increases rapidly. Emis- 
sion from these regions is, however, negligible in the XRT 
band due to the relatively low temperatures there. After shock 
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Fig. 11. — The logarithmic temperature in Kelvin for model m2rlcold at three post-breakout times: 68.2 s, 124.9 s, 202.3 s, from left to right. The white 
contour shows the location of the thermalization depth at 0.75 keV. The line-of-sight is along the equator. Prior to shock eruption, the thermalization depth is 
located in the transition layer between the wind and the progenitor model. Post-shock-breakout, the thermalization depth is in the shocked wind, ahead of the 
reverse shock. 
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Fig. 12.— XRT band (0.1 - 10 keV) light curves for model m2rlcold at 
viewing angles of 0, 45, and 90 degrees. The top panel shows the X-ray lu- 
minosity light curves and the bottom panel shows the XRT count rate light 
curves corrected for detector response and absorption along the line of sight. 
The observed XRT light curve of XRO 080109 is also shown for comparison. 
The light curves are double-peaked in shape, with the first peak correspond- 
ing to emission from the spot where the shock is erupting from the progenitor 
surface and the second peak corresponding to enhanced emission from the 
equator as the shocks cross. The second peak, which has a much larger lu- 
minosity than the first, is characterized by softer emission relative to the first 
and so results in a similar X-ray count rate to the first peak, once corrections 
for absorption are made. The light curve for a viewing angle of 90 degrees 
(along the axis of symmetry) is generally less luminous than the other view- 
ing angles, even early on. This is because once the shocks erupt from the 
poles, they expand quickly and obscure the very bright ring on the surface of 
the progenitor where the shock continues to erupt. 

breakout the thermalization layer lies in between the forward 
and reverse shocks, the region comprise of shocked, acceler- 
ated wind material. 

Light curves of explosion models m2rlcold and m2rlhot 
(the smaller progenitor) are shown in Figures [12] and [13] for 
three viewing angles: along the equator (0 degrees), 45 de- 
grees, and along the axis of symmetry (90 degrees). The 
widths of the light curves are roughly 100 seconds at all an- 
gles for both models and the total X-ray burst time is around 
500 to 600 seconds. These values are very close to what was 
observed for SN 2008D, shown in Figure [8] Examining the 
shape of the model light curves shows that for degree view- 
ing angles the light curves rise very quickly. This is due to two 
effects. The first is simply because the emission from both 



Fig. 13.— XRT band (0.1 - 10 keV) light curves for m2rlhot at viewing 
angles of 0, 45, and 90 degrees. The top panel shows the X-ray luminosity as 
a function of time and the bottom panel shows the predicted XRT count rate 
accounting for detector response and absorption due to neutral matter. The 
observed XRT light curve of XRO 080109 is also shown for comparison. 
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Fig. 14. — XRT band (0.1 - 10 keV) light curves for m7r6cold at viewing 
angles of 0, 45, and 90 degrees. The top panel shows the X-ray luminosity as 
a function of time and the bottom panel shows the predicted XRT count rate 
accounting for detector response and absorption due to neutral matter. The 
observed XRT light curve of XRO 080109 is also shown for comparison. 
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Fig. 15. — XRT band (0.1 - 10 keV) light curves for m7r6hot at viewing 
angles of 0, 45, and 90 degrees. The top panel shows the X-ray luminosity as 
a function of time and the bottom panel shows the predicted XRT count rate 
accounting for detector response and absorption due to neutral matter. The 
observed XRT light curve of XRO 080109 is also shown for comparison. 



poles is visible. The second is due to the nature of an aspheri- 
cal breakout and the effects of light travel corrections. As the 
bipolar shocks continue to erupt from the progenitor surface, 
the brightest emission is coming from where the shocks are 
just reaching the surface (essentially two rings moving across 
the stellar surface from the poles to the equator). Thus, the 
brightest emitting regions are moving rapidly toward the ob- 
server causing a pile up of emission once light travel time 
corrections are made. This effect is reduced, or eliminated, at 
higher viewing angles because the brightly emitting rings are 
no longer moving so much toward the observer. 

Table |2] lists the resulting light curve characteristics for our 
models m2rlcold and m2rlhot, as well as m2rlsph for com- 
parison. As expected, the spherical explosion has a higher 
peak luminosity, but a much shorter FWHM and overall burst 
time, A^ The width of the light curve for m2rlsph is set by 
the light crossing time of the progenitor. The simulated XRT 
count rates for models m2rlhot and, especially, m2rlcold are 
very similar to those we find for XRO 080109/SN 2008D. 
The light curve rise time and FWHM for model m2rlcold at 
a viewing angle of 45 degrees (red dashed curve in Figure [12]) 
is a good match to SN 2008D. 

The X-ray spectra for models m2rlcold and m2rlhot are 
shown in Figure [17] These spectra are corrected for the de- 
tector response function and for X-ray absorpotion (assuming 
A^H = 1-7 X 10^^). They are also averaged over the burst time, 
A^ as was done for t he observations of XRO 080109/SN 
2008D (see Figure |9| .Soderberg et al. 
2009| ). The shapes of the spectra at t 
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le different viewing 

Indeed there is not much differ- 

The time-averaging 



angles are very similar. 

ence between the two different models 

washes away the major differences that are apparent in the 

light curves. The shape of the spherical explosion spectra are 

also very similar. They are generally brighter, but this is due 

to a shorter averaging time. The spectra in each case are softer 

than the spectrum of XRO 080109 (see Figure[9]). 

The light curves for the simulations in the larger progenitor 
(Figures 14 and [TSj are characterized by much longer time 
scales and overall less bright emission. The FWHM for mod- 
els m7r6cold and m7r6hot range from around 100 to 200 sec- 
onds at various angles while the total burst times are from 300 
to 1000 seconds (see Table[2]). The peak luminosities and total 
radiated energies are less than in the analogous explosions in 
the smaller progenitor. While the ratio of explosion energies 
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Fig. 16. — Instantaneous XRT spectra for model m2rlcold at four different 
times during the first 150 seconds of breakout emission. The times of the 
spectra are indicated in the legend. The spectra show significant evolution to 
softer emission as the burst proceeds. 
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Fig. 17. — Time-integrated X-ray spectra for the smaller progenitor ex- 
plosion models, m2rlcold (top) and m2rlhot (bottom) at different view- 
ing angles. The spectra are corrected for detector response and absorp- 
tion along the line of sight, assuming a neutral matter column depth of 
A^H = 1-7 X 10^^ cm~^. For comparison, we also plot the spectrum for the 
spherical explosion, m2rlsph, and the observed XRT spectra of XRO 080109. 

to ejecta masses are roughly equivalent across all simulations, 
the reduced luminosities in the larger progenitor can be ex- 
plained by a slightly lower shock velocity during the burst 
and a lower wind density at the radius of shock breakout (i.e., 
the radius of t he p rogenitor star). As we discuss in more de- 
tail in Section [421 the density of the wind plays an important 
role in the strength of the X-ray emission because the thermal- 
ization depths during the bursts lie in the shocked wind. The 
light curve for m7r6cold is similar in shape to m2rlcold and 
m2rlhot but with longer timescales. Model m7r6hot, how- 
ever, exhibits a dramatically double-peaked light curve. 

As discussed in Section [2] we ran simulations of model 
m2rlcold at three different resolutions. Figure [19] shows the 
X-ray light curves for these three simulations plotted together. 
The light curves of the low- and fiducial resolution simula- 
tions are extremely similar. The high-resolution case varies 
only slightly from the lower-resolution light curves. The early 
peak is brighter, du e to the slightly greater shock speeds just 
after breakout (see |4.2| ). The second peak is slightly dim- 
mer, due to the lesser degree of extension of the southern 
shock structure in the high-resolution simulation, leading to 
a smaller emitting area as seen from a viewing angle of de- 
grees. 

4.1. Dependence on Metallicity and X-ray Absorption 

The absorptive opacities depend strongly on the metallic- 
ity of the absorbing gas. To illustrate this we have calculated 
simulated spectra and light curves of our explosion models 
using metallicity values of O.IZ© and O.OZ© (i.e., metal-free). 
Figures [20] and [2T] show the light curves for model m2rlcold 
for metallicities of Z = 0. IZ© and Z = 0. The spectra for these 
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Fig. 18. — Time-integrated X-ray spectra for the larger progenitor ex- 
plosion models, m7r6cold (top) and m7r6hot (bottom) at different view- 
ing angles. The spectra are corrected for detector response and absorp- 
tion along the line of sight, assuming a neutral matter column depth of 
Nil = l.l X 10^^ cm~^. For comparison, we also plot the spectrum for the 
spherical explosion, m7r6sph, and the observed XRT spectra of XRO 080109. 



Fig. 20.— XRT band (0.1 - 10 keV) light curves for m2rlcold at view- 
ing angles of 0, 45, and 90 degrees calculated assuming a metallicity of 
Z = 0. 1Z0 . The top panel shows the X-ray luminosity as a function of time 
and the bottom panel shows the predicted XRT count rate accounting for de- 
tector response and absorption due to neutral matter. The observed XRT light 
curve of XRO 080109 is also shown for comparison. 
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Fig. 19. — XRT count rate light curves for model m2rlcold at three differ- 
ent resolutions. 



cases are shown in Figure 22 The lower metalHcity drives the 
thermaUzation depth deeper into the explosion where the tem- 
peratures are higher, resulting in significantly increased emis- 
sion (note the difference in scale in Figure [21]). For the case 
of metal-free gas, this increase is dramatic. Due to the signifi- 
cantly reduced absorptive opacities, the thermalization depths 
are pushed down below the reverse shock and into the deep, 
very hot regions of the ejecta. The calculated emission for 
this metallicity is orders of magnitude greater than the other 
models and the observations of SN 2008D. The difference in 
the emission characteristics between our fiducial models with 
Z = O.SZq and the models with Z = 0. IZq is less drastic. This 
is because for Z = O.IZ©, the thermalization depths are still 
above the reverse shock and the temperature in between the 
forward and reverse shocks does not vary greatly. The metal- 
free case spectra also noticeably lack the deep 'absorption' 
features at ~ 0.4 keV and ~ 0.9 keV. These 'absorption' fea- 
tures are caused by a significant increase in the absorptive 
opacities at these energies (due to the presence of metals), 
pushing the thermalization depths to larger, cooler radii. We 
note that our fiducial metallicity of Z = O.SZq is consistent, 
within the accuracies, with three independent measurements 
of the metal l icity of the region a r ound SN 2008D: [Soderberg 
eraL| ( [2Q08l ); |Thone et al. | ( [20091 ) ; |Modjaz ^t al.| ( |2010| ). 

We plot XRT count rates and spectra for model m2rlcold 
using different values of A^h in Figures [23] and [24) If the loca- 
tion of SN 2008D in its host galaxy were particularly dense. 
X-ray absorption by neutral matter atoms in the host galaxy 
may be significant, warranting a A^h beyond the Galactic value 
of 1.7x10^^ cm~^. Increased hydrogen column depth in- 
creases the absorption of lower-energy X-ray photons. Since 
our simulated spectra for the non-zero metallicity cases are 
dominated by emission from below about 1 keV, increasing 




Fig. 21.— XRT band (0.1 - 10 keV) light curves for m2rlcold at viewing 
angles of 0, 45, and 90 degrees calculated assuming a metallicity of Z = OZq , 
metal-free gas. The top panel shows the X-ray luminosity as a function of 
time and the bottom panel shows the predicted XRT count rate accounting 
for detector response and absorption due to neutral matter. 

A^H dramatically reduces the resultant count rates. Due to 
a significantly greater amount of hard emission in the Z = 
model, the light curve and spectrum are effected very little by 
an increase in A^h- It is possible to increase the neutral matter 
column depth so that the peak count rates for the metal-free 
model are similar to those observed for XRO 080109 (< 10 
cnts s~^). This requires, however, an unrealistically high value 
of A^H, greater than 10^^ cm~^ and results in a very narrow light 
curve as the softer emission after the first 20 seconds is almost 
entirely absorbed. 

4.2. Color Temperature Enhancement 

In our simulations, the thermalization depths where the X- 
ray spectrum is formed, lie in the region of shocked wind, be- 
hind the forward shock and ahead of the reverse shock. Since 
we assume local emission with a black-body spectrum, the re- 
sultant luminosity is a strong function of the temperature used 
in calculating the spectrum. In the region in between the for- 
ward and reverse shock, this temperature is dependent on the 
forward shock velocity and the wind density. To see this, as- 
sume that enthalpy flux is conserved at the shock front such 
that, in the frame of the shock, 

PiVi{\v\ + ei+Pi/ pi) = p2V2C2vl + e2+P2/ pi), (7) 
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Fig. 22. — Time-integrated X-ray spectrum for model m2rlcold using 
Z = 0. 1Z0 (top) and Z = (bottom) along with the spectrum for the spherical 
explosion at the respective metallicities. The observed XRT spectrum of XRO 
080109 is also plotted. 
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Fig. 23. — XRT count rate light curve (top) and time-integrated XRT 
spectrum (bottom) for m2rlcold calculated using metallicities of 0.1 Z© and 
O.5Z0 and neutral matter column depth of 2 x 10^^ cm~^ and 5 x 10^^ cm~^. 
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Fig. 24. — XRT count rate light curve (top) and time-integrated XRT spec- 
trum (bottom) for m2rlcold calculated for metal-free gas and neutral matter 
column depth of 2 x 10^^ cm~^ and 5 x 10^^ cm~^. 

where p is gas density, v is gas velocity, e is the specific inter- 
nal energy of the gas, P is the gas pressure, and subscript 1 de- 
notes pre-shock values and subscript 2, post-shock values. We 
can assume that the post-shock internal energy is dominated 
by contributions from radiation, 62 = '^Pi/ pi = <3rad^2^/P2, and 
the density jump at the shock is p2 = 1 pi, for a strong shock. 
So long as the shock is strong, the upstream internal energy 
is negligible, ei ^ Pi ^0. Conservation of mass flux at the 



shock also gives V2 = vipi/ p2 = vi/7, in the shock frame. In 
the frame of the progenitor star, vi is the shock speed, Vsh, as 
long as the shock is moving fast relative to the pre-shock gas. 
Substituting these relations into equation ^ and solving for 
the post-shock temperature yields 

T2^1.6xlO'p[{\vl{^,,K, (8) 

g cm~^ and 
^^^ cm s~^ appropri- 



where p-n is the wind density in units of 10 

Vsh, 10 is the shock velocity in units of 10 

ate for shock breakout into a Wolf-Ray et wind. Thus, the 

post- shock temperature may be increased by enhancing the 

wind density or the energy of the explosion (which, in turn, 

increases Vsh). 

In order to demonstrate the influence that an increased post- 
shock temperature has on our simulated spectra and light 
curves, we present the spectrum and light curve for model 
m2rlcold calculated by using a temperature that had been en- 
hanced by a factor of 1.8 above the temperature found in our 
hydrodynamic simulations. Such an increase in tem perature 



would result from a factor of ten increase in piv^^. [Tanaka 
|et al.| p009) find that one dimensional explosion models with 
kinetic energy to ejecta mass ratios of 1.4 - 1.7 (10^^ erg/M©) 
fit the late time spectrum and light curve of SN 2008D best. 
This ratio is about a factor of 2 greater than what we have 
used in our most simulations, including m2rlcold. This alone 
could account for a ten-fold increase in the shock ram pressure 
piVgjj because the mass averaged- velocity will scale roughly 
as Vgj (X E]^/Mqj and the shock speed will far exceed Vej- 

The results of this enhanced temperature calculation for 
model m2rlcold are shown in Figure 25 The simulated light 
curves and spectra are all for a metallicity of O.SZq and a 
viewing angle of degrees. Figure 25 shows the behav- 
ior of the simulated emission with increased A^h- The en- 
hanced temperature calculation yields a peak X-ray luminos- 
ity of 4.9 X 10"^^ erg s~^ and radiates a total of 2.0 x 10"^^ ergs. 
For the case of A^h = 1-7 x 10^^ cm~^, the X-ray count rate 
peaks at 65 s~\ much higher than was observed for XRO 
080109. Increasing the column depth brings down the peak 
count rate, narrows the width of the light curve, and hardens 
the spectrum. For A^h = 2 x 10^^ cm~^, the peak count rate 
is 10 s~\ slightly greater than the observed value. The light 
curve at this column depth is somewhat narrower and rises 
more quickly than that of XRO 080109. The spectrum is also 
still too soft and devoid of significant count rates above 2 keV. 
In Figure [26] we show the angular dependence of the sim- 
ulated light curve and spectra for the enhanced temperature 
version of m2rlcold calculated with a metallicity of O.SZq 
and A^H = 1-7 X 10^^ cm~^. The spectrum is harder than the 
fiducial models, computed without temperature enhancement, 
but is still softer than that of XRO 080109. 

The color temperature of the emission may also be en- 
hanced by non-LTE effects. The radiation escaping from the 
shock during breakout from a Wolf-Rayet star is not in ther- 
mal equilibrium wit h the gas at r* = 2/3 ( |Katz et aL]|2009 
Nakar & Sari||201Q| ). In fact, during the early parts of the 



shock breakout the radiation temperature may be or ders of 
magnitude greater than the gas temperature at r* = 2/3 (Nakar 
,& Sari,2010j . Over the course of the first tens of seconds of 
the X-ray burst, the radiation temperature will quickly drop 
down closer to the gas temperature. This time-dependent en- 
hancement of the radiation temperature, if accounted for in 
our calculations, could produce the hard X-rays that are miss- 
ing in our spectra while not dramatically increasing the lumi- 
nosity at later times in the light curve. 
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Fig. 25. — XRT count rate light curves (top) and time-integrated XRT 
spectra (bottom) for m2rlcold calculated using a temperature enhancement 
factor of 1.8 at various values of A^h- 
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Fig. 26. — XRT count rate light curve (top) and time-integrated XRT spec- 
tra (bottom) for the enhanced temperature version of m2rlcold at various 
viewing angles assuming a neutral matter column depth of 2 x 10^^ cm~^. 
Also shown are the observed XRT light curve (top) and spectrum (bottom). 

4.3. SED Shape 

The shapes of the spectral energy distributions of our sim- 
ulated shock breakout models are not well-fit by single- tem- 
perature and radius spherical black-bodies. This is the case 
for both the jet-driven explosion simulations and the spherical 
explosion models. Figures 27 and 28 show the X-ray spectra 
of models m2rlcold and m2rlsph along with spherical single- 
temperature and radius black-body SEDs corrected for XRT 
detector response and X-ray absorption. As is shown, no sin- 
gle temperature and radius black-body can fit both the soft and 
hard parts of the spectrum in either the jet-driven or spherical 
explosions. This demonstrates that, regardless of shock break 
out geometry, such simple black-body models are not applica- 
ble. Shock breakout is too dynamic a process to be modeled 
by a single average radius and a single average temperature 
black-body, even for spherical breakouts. Simple dynamical 
models with radii expanding at constant velocities and tem- 
peratures cooling adiabatically could provide vast improve- 
ments over simple static models and could provide a better 
relation to the physical process of shock breakout. 

5. DISCUSSION AND CONCLUSIONS 

We have modeled the emission arising from the breakout 
of an aspherical supernova shock assuming that the emission 
is thermal in nature. Our hydrodynamic simulations do not 
account for radiative effects such as disparate ion and radia- 
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Fig. 27. — X-ray spectrum of jet explosion model m2rlcold (black line) 
plotted with two spherical, single temperature and radius black-body SEDs, 
corrected for XRT detector response and X-ray absorption. The red dash- 
dot line is a black-body SED with spherical radius of 1.2 x 10^^ cm and 
temperature 3 x 10^ K. The blue dashed line is a black-body SED with radius 
of 8 X 10^^ cm and temperature 7 x 10^ K. No single spherical black-body 
SED fits the entire spectrum of m2rlcold well. 
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Fig. 28. — X-ray spectrum of spherical explosion model m2rlsph (black 
line) plotted with two spherical, single temperature and radius black-body 
SEDs. The red dash-dot line is a black-body SED with a spherical radius of 
5 X 10^^ cm and temperature 5 x 10^ K. The blue dashed line is a black-body 
SED with a radius of 7 x 10^^ cm and temperature 1 x 10^ K. Despite being 
a spherical explosion, no single spherical black-body SED fits the spectrum 
well. This is because the radius of the emitting region in the explosion model 
is changing rapidly and no single color temperature can describe the emission 
accurately. 



tion temperatures, radiative pre- acceleration of gas ahead of 
the shock, escape of photons from the shock at optical depths 
greater than 1, etc. Despite the limitations of our simulations, 
we have demonstrated that in the more general case of a non- 
spherical explosion, the breakout emission is dramatically dif- 
ferent from that expected from a spherically- symmetric shock 
breakout. One of the most important results of this work is 
that the timescales of the light curves for aspherical shock 
breakout are not set by the light travel time across the pro- 
genitor star, but are instead related to the much longer shock 
crossing time of the progenitor. Thus, shock breakout light 
curves contain information about the geometry of the shock 
structure as well as the radius of breakout. We also show that 
for aspherical shock breakouts, the observer's viewing angle 
can play an important role in determining the shape of the 
observed light curve. These general results apply to any, ar- 
bitrarily aspherical supernova, not just a jet-driven, bipolar 
supernova. 

The fiducial emission models are generally under-luminous 
and have soft spectra compared to the actual observations. 
Tanaka et al.| ([2009 ) advocate explosions with higher ener- 
gies than what we have simulated and S oderberg et al.| ( [2008] ) 
posit that the wind around the progenitor of SN 2008D was 
optically-thick, thus very dense. Either an increased explo- 
sion energy or increased wind density could increase the X- 
ray luminosity of shock breakout; however, even if the wind 
and explosion energies were increased to give luminosities 
commensurate with those inferred for XRO 080109, the spec- 
trum would still remain lacking in sufficient luminosity above 
2 keV to match the observed spectrum (see Figure [25]). This 
seems to indicate the need for a scattering of the therma l X- 
ray photons to hi gher energie s , as p rescribed by ,Soderberg| 
et al.| p008| ) and |Wang et al.] ( |2008| ), or other, non-thermal 
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jaz et al.| ( |2QQ8| ) for several other Type Ib/c supemovae with 
double-peaked oxygen lines. ^ A double-peaked oxygen line 
may also result from an asymmetric distribution of radioac- 
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Fig. 29. — Nickel (left) and oxygen (right) mass fractions for model 
m2rlcold at the end of the simulation (10^ seconds). Both distributions show 
extension along the symmetry axis with a great deal of small scale, turbulent 
structure. The oxygen distribution also shows cavities caused by the jets. 



processes (Katz et al. 2009: Nakar & SaripOlOl ). This could 
serve to harden our simulated spectra and also lengthen the 
burst time since, in our models, the later parts of the light 
curve are made up of a larger fraction of soft, highly-absorbed 
photons. 

The simulated spectra may also be hardened by relaxing 
the assumption of thermal equilibrium between the matter 
and radiation. During shock breakout from a massive, com- 
pact Wo lf-Rayet pr ogenitor the matte r and radiation are not 
in LTE dKatz et al . 2009; N akar & S ari 2010). Our method 
does no t account for non-LTE effects. As discussed by |Nakar| 
|& Sari| (|2010), if the matter and radiation are not in LTE, 
the observed color temperature could be much greater than 
the matter temperature at r* ^ 2/3. Because the internal en- 
ergy would not be similarly enhanced above the LTE case, 
the bolometric luminosities for non-LTE breakouts will not 
depend strongly on the coupling between matter and radia- 
tion, as the color temperature does. The color- temperature 



enhancing eff ects in non-LTE breakouts discussed by |Nakar 
& Sari| ( |2010| ) are strongly time-dependent with the greatest 



difference in the temperatures of radiation and matter being 
at the instant o f shock brea kout. If non-LTE effects, such 
as described by |Nakar & Sa ri (2010), were included in our 
simulations, the spectra would be significantly hardened, es- 
pecially at early times in the breakout, and the X-ray lumi- 
nosities modestly enhanced. This could account for the miss- 
ing hard X-ray emission in our simulated spectra and increase 
the XRT count rates without the need to invoke a very dense 
wind or a greatly enhanced exp losion energy. We note, how- 
ever, that [Nakar & Sari| ( |2010| ) do not include the effects of 
a wind surrounding the progenitor. In our calculations, the 
presence and character of the wind is an integral factor in the 
breakout emission formation since the thermalization depths 
in our simulations lie at, or ahead of, the contact discontinuity 
between the ejecta and the wind. 

The presence of a double-peaked oxygen line in the spec- 
tra of SN 2008D at about 109 days may be evidence for as- 
phericity in the explosion (Modj az et al.|2009) . [Modjaz et al.| 
([2009 ) suggest that this may be evidence for a ring-like distri 
bution of oxygen, similar to the conclusions reached by|Mod 



tive cobalt powering the excitation of oxygen ( Modjaz et al. 
"2008, Gera rdy et al.|2 000). If the cobalt distribution were as- 
pherical, however, other radioactively-excited lines would be 
expected to show a similar double-peaked behavior, which is 
not the case for SN 2 08D or the Typ e Ib/c SNe discussed in 
[Modjaz et aT] ( [2008] ). [Tanaka et aL] (|2009 ) find that a frac- 
tion of nickel must be mixed outward to adequately model 
the spectra of SN 2008D, further indicating asphericity of the 
explosion. Figure 29 shows the nickel and oxygen mass frac- 
tions from our model m2rlcold at 10^ seconds. In our simula- 
tions, both the nickel and oxygen distributions are aspherical, 
as is shown in Figure [29] This may be able to account for 
the double-peaked oxygen line in the spectra of SN 2008D, 
however detailed spectral synthesis calculations are needed 
to be sure. Not e that our oxygen distribution is prolate, not 
ring-like as Mo djaz et al. ( 2008| ) recommend for explaining 
double-peaked oxygen lines in SNe. 

The distribution of intermediate-mass elements such as 
oxygen we find in our simulations differs from what some 
other g roups find in similar studies of aspherical SNe. |Maeda[ 
et al.[ ( [2002[ ) present hydrodynamic simulations of aspher- 



ical CCSNe targeted to explaining the observations of SN 
1998bw. They show that the intermediate elements are ejected 
from the explosion in an equatorial torus. We find that the in- 
termediate elements are ejecta in a bipolar geometry, as shown 
in Figure [29] The difference in the two results comes from 
the manner in which the explosions are initiated. In |Maeda[ 
et al. , the explosions are started by depositing kinetic energy 
asymmetrically in the center of the progenitor. This pushes 
the intermediate elements outward while simultaneously the 
more energetic material near the poles compresses them into a 
toroidal geometry. In our simulation the explosions are driven 
entirely by the bipolar jets. The intermediate-mass material 
near the progenitor's equator is accreted into the central en- 
gine. Intermediate elements are entrained in the jets and car- 
ried out into a configuration that resembles that of the jets 
themselves. The equatorial torus in our simulations is com- 
prised primarily of helium. The final distribution of interme- 
diate mass elements thus depends on the mode of asymmet- 
ric energy input. Determining the distribution observationally 
may thus help to cons train models. 

2009) present early spectropolarimetric ob- 



Maund et aL[ 

servations of SN 2008D. They find that the continuum polar 
ization is relatively small, indicating that the supernova pho- 
tosphere may be only slightly aspherical. They also find that 
there is significant polarization in certain spectral lines, indi- 
cating that the line-forming regions of various elements are 
markedly aspherical. This is in qualitative agreement with 
the results of our simulations. The photospheres at late times 
are nearly round, while the detailed composition structure is 
dramatically asymmetric, as shown in Figure 29 A late-time 
shock structure that is nearly round, or at least not dramat- 
ically as pherical, is also consistent with the radio measure- 
ments of Bietenholz et al. (2009). 

The number of observed supernova shock breakouts has 
been increasing, and this trend is likely to continue, opening 
a new window for exploring core-collapse supernova through 

^ This interpretation of double-peaked oxygen lines is somewhat contro- 
versial. Milisavljevic et al. (2010) argue that this feature is caused by an 
oxygen doublet at these wavelengths. 
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observation. A multidimensional interpretation of these ob- 
servations will be critically important to gaining an accu- 
rate understanding of supernova shock breakout. We find 
that aspherical shock breakout can account for the light curve 
time scales observed for XRO 080109/SN 2008D without re- 
quiring an extremely dense wind or abnormally large Wolf- 
Rayet progenitor star. Our models with the smaller progen- 
itor, m2rlcold and m2rlhot, have light curves with FWHMs 
of about 100 seconds, roughly matching XRO 080109. We 
find that our models are generally under-luminous and have 
spectra too soft to match XRO 080109; however this is likely 
due to our assumption of LTE. The shock breakout from a 
WR progenitor is not in LTE and the radiation temperature in 
the non-LTE case c ould be significan tly enhanced above the 
matter temperature (Nakar & Sari 2010 ). Accounting for non- 
LTE effects, then, could brighten our simulated light curves 
and harden the simulated spectra. Inclusion of non-LTE ef- 
fects will be presented in future work. 



The authors thank Peter Hoflich and Ehud Nakar for help- 
ful discussion concerning the physics of shock breakout. 
S.M.C. is very grateful for support provided through NASA 
Earth and Space Science Fellowship program. J.C.W. is sup- 
ported in part by NSF grant AST-0707769 and in part by 
DOE funding to the Center for Laser Experimental Research 
(CLEAR) through the University of Michigan. The soft- 
ware used in this work was in part developed by the DOE- 
supported ASC/Alliance Center for Astrophysical Thermonu- 
clear Flashes at the University of Chicago. The authors ac- 
knowledge the Texas Advanced Computing Center (TACC) 
at The University of Texas at Austin for providing high- 
performance computing, visualization, and data storage re- 
sources that have contributed to the research results reported 
within this paper. 



REFERENCES 



Almgren, A. S., Bell, J. B., Rendleman, C. A., & Zingale, M. 2006, ApJ, 637, 

922 
Berger, E. & Soderberg, A. M. 2008, GRB Coordinates Network, 7159, 1 
Bietenholz, M. R, Soderberg, A. M., & Bartel, N. 2009, ApJ, 694, L6 
Blondin, J. M., Mezzacappa, A., & DeMarino, C. 2003, ApJ, 584, 971 
Buras, R., Janka, H., Rampp, M., & Kifonidis, K. 2006, A&A, 457, 281 
Burrows, A., Dessart, L., Livne, E., Ott, C. D., & Murphy, J. 2007, ApJ, 664, 

416 
Burrows, A., Livne, E., Dessart, L., Ott, C. D., & Murphy, J. 2006, ApJ, 640, 

878 
Chevalier, R. A. & Fransson, C. 2008, ApJ, 683, L135 
Colgate, S. A. 1968, Canadian Journal of Physics, 46, 476 
— . 1974, ApJ, 187,333 

Couch, S. M., Wheeler, J. C, & Milosavljevic, M. 2009, ApJ, 696, 953 
Dickey, J. M. & Lockman, F. J. 1990, ARA&A, 28, 215 
Ensman, L. & Burrows, A. 1992, ApJ, 393, 742 
Fryxell, B., Arnett, D., & Mueller, E. 1991, ApJ, 367, 619 
Fryxell, B., Olson, K., Ricker, R, Timmes, R X., Zingale, M., Lamb, D. Q., 

MacNeice, R, Rosner, R., Truran, J. W., & Tufo, H. 2000, ApJS, 131, 273 
Gerardy, C. L., Pesen, R. A., Hoflich, R, & Wheeler, J. C. 2000, AJ, 119, 

2968 
Gorosabel, J., de Ugarte Postigo, A., Castro-Tirado, A. J., Agudo, L, Jelinek, 

M., Leon, S., Augusteijn, T., Fynbo, J. P. U., Hjorth, J., Michalowski, 

M. J., Xu, D., Ferrero, R, Kann, D. A., Klose, S., Rossi, A., Madrid, J. R, 

LLorente, A., Bremer, M., & Winters, J. 2008, ArXiv e-prints 
Hoflich, P., Khokhlov, A., & Wang, L. 2001, in American Institute of 

Physics Conference Series, Vol. 586, 20th Texas Symposium on relativistic 

astrophysics, ed. J. C. Wheeler & H. Martel, 459-+ 
Katz, B., Budnik, R., & Waxman, E. 2009, ArXiv e-prints 
Khokhlov, A. M., Hoflich, R A., Oran, E. S., Wheeler, J. C, Wang, L., & 

Chtchelkanova, A. Y. 1999, ApJ, 524, L107 
Kifonidis, K., Plewa, T., Janka, H.-T., & MuUer, E. 2003, Astronomy & 

Astrophysics, 408, 621 
Kifonidis, K., Plewa, T., Scheck, L., Janka, H.-T., & MuUer, E. 2006, 

Astronomy & Astrophysics, 453, 661 
Klein, R. I. & Chevalier, R. A. 1978, ApJ, 223, L109 
Kong, A. K. H. & Maccarone, T. J. 2008, The Astronomer's Telegram, 1355, 

1 
Liska, R. & Wendroff, B. 2003, SIAM J. Sci. Comput., 25, 995 
Maeda, K., Nakamura, T., Nomoto, K., Mazzali, P. A., Patat, F, & Hachisu, 

I. 2002, ApJ, 565, 405 
Magee, N. H., Abdallah, Jr., J., Clark, R. E. H., Cohen, J. S., Collins, L. A., 

Csanak, G., Fontes, C. J., Gauger, A., Keady, J. J., Kilcrease, D. P., & 

Merts, A. L. 1995, in Astronomical Society of the Pacific Conference 

Series, Vol. 78, Astrophysical Applications of Powerful New Databases, 

ed. S. J. Adelman & W. L. Wiese, 51-h 
Matzner, C. D. & McKee, C. R 1999, ApJ, 510, 379 
Maund, J. R., Wheeler, J. C, Baade, D., Patat, R, Hoflich, R, Wang, L., & 

Clocchiatti, A. 2009, ApJ, 705, 1139 



Mazzali, P. A., Valenti, S., DeUa Vafle, M., Chincarini, G., Sauer, D. N., 
Benetti, S., Plan, E., Piran, T., D'Elia, V, EHas-Rosa, N., Margutti, 
R., Pasotti, R, Antonelli, L. A., Bufano, R, Campana, S., Cappellaro, 
E., Covino, S., D'Avanzo, P., Piore, P., Pugazza, D., Gilmozzi, R., 
Hunter, D., Maguire, K., Maiorano, E., Marziani, P., Masetti, N., Mirabel, 
P., Navasardyan, H., Nomoto, K., Palazzi, E., Pastorello, A., Panagia, 
N., Pellizza, L. J., Sari, R., Smartt, S., Tagliaferri, G., Tanaka, M., 
Taubenberger, S., Tominaga, N., Trundle, C, & Turatto, M. 2008, Science, 
321,1185 

Milisavljevic, D., Pesen, R. A., Gerardy, C. L., Kirshner, R. P., & Challis, P. 
2010, ApJ, 709, 1343 

Modjaz, M., Bloom, J. S., Pilippenko, A. V, Kewley, L., Perley, D., & 
Silverman, J. M. 2010, ArXiv e-prints 

Modjaz, M., Kirshner, R. R, Blondin, S., Challis, R, & Matheson, T. 2008, 
ArXiv e-prints 

Modjaz, M., Li, W., Butler, N., Chornock, R., Perley, D., Blondin, S., 
Bloom, J. S., Pilippenko, A. V, Kirshner, R. P., Kocevski, D., Poznanski, 
D., Hicken, M., Foley, R. J., Stringfellow, G. S., Berlind, R, Barrado y 
Navascues, D., Blake, C. H., Bouy, H., Brown, W. R., ChaUis, R, Chen, 
H., de Vries, W. H., Dufour, P., Falco, E., Friedman, A., Ganeshalingam, 
M., Garnavich, P., Holden, B., Illingworth, G., Lee, N., Liebert, J., Marion, 
G. H., Olivier, S. S., Prochaska, J. X., Silverman, J. M., Smith, N., Starr, 
D., Steele, T. N., Stockton, A., WiUiams, G. G., & Wood-Vasey, W. M. 
2009, ApJ, 702, 226 

Nakar, E. & Sari, R. 2010, ArXiv e-prints 

Obergaulinger, M., Aloy, M. A., & Muller, E. 2006, A&A, 450, 1 107 

Plewa, T. & Muller, E. 1999, A&A, 342, 179 

Rybicki, G. B. & Lightman, A. P. 1986, Radiative Processes in Astrophysics 
(Wiley-VCH) 

Soderberg, A. M., Berger, E., Page, K. L., Schady, P., Parrent, J., Pooley, D., 
Wang, X., Ofek, E. O., Cucchiara, A., Rau, A., Waxman, E., Simon, J. D., 
Bock, D., Milne, P. A., Page, M. J., Barentine, J. C, Barthelmy, S. D., 
Beardmore, A. P., Bietenholz, M. R, Brown, P., Burrows, A., Burrows, 
D. N., Byrngelson, G., Cenko, S. B., Chandra, P., Cummings, J. R., Fox, 
D. B., Gal- Yam, A., Gehrels, N., Immler, S., Kasliwal, M., Kong, A. K. H., 
Krimm, H. A., Kulkarni, S. R., Maccarone, T. J., Meszaros, P., Nakar, E., 
O'Brien, P. T., Overzier, R. A., de Pasquale, M., Racusin, J., Rea, N., & 
York, D. G. 2008, Nature, 453, 469 

Suzuki, A. & Shigeyama, T. 2010, ApJ, 717, L154 

Tanaka, M., Tominaga, N., Nomoto, K., Valenti, S., Sahu, D. K., Minezaki, 
T., Yoshii, Y, Yoshida, M., Anupama, G. C, Benetti, S., Chincarini, G., 
Delia Valle, M., Mazzali, R A., & Plan, E. 2009, ApJ, 692, 1131 

Thone, C. C, Michalowski, M. J., Leloudas, G., Cox, N. L. J., Fynbo, J. R U., 
SoUerman, J., Hjorth, J., & Vreeswijk, R M. 2009, ApJ, 698, 1307 

Wang, L., Howell, D. A., Hoflich, R, & Wheeler, J. C. 2001, ApJ, 550, 1030 

Wang, L. & Wheeler, J. C. 2008, ARA&A, 46, 433 

Wang, X., Li, Z., Waxman, E., & Meszaros, P. 2008, in American Institute 
of Physics Conference Series, Vol. 1065, American Institute of Physics 
Conference Series, ed. Y.-R Huang, Z.-G. Dai, & B. Zhang, 279-284 

Wheeler, J. C, Meier, D. L., & Wilson, J. R. 2002, ApJ, 568, 807 

Wheeler, J. C, Yi, I., Hoflich, R, & Wang, L. 2000, ApJ, 537, 810 

Woosley, S. E., Langer, N., & Weaver, T. A. 1995, ApJ, 448, 315 



1 6 ASPHERIC AL SHOCK BREAKOUT 

Zingale, M., Dursi, L. J., ZuHone, J., Calder, A. C, Fryxell, B., Plewa, T., 
Truran, J. W., Caceres, A., Olson, K., Ricker, P. M., Riley, K., Rosner, R., 
Siegel, A., Timmes, F. X., & Vladimirova, N. 2002, ApJS, 143, 539 



